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ABSTRACT 

The RAPid Telescopes for Optical Response (RAPTOR) system at Los Alamos National Laboratory 
observed GRB 050319 starting 25.4 seconds after 7-ray emission triggered the Burst Alert Telescope 
(BAT) on-board the Swift satellite. Our well sampled light curve of the early optical afterglow is 
composed of 32 points (derived from 70 exposures) that measure the flux decay during the first hour 
after the GRB. The GRB 050319 light curve measured by RAPTOR can be described as a relatively 
gradual flux decline (power-law index a = —0.37) with a transition, at about ~ 400 s after the GRB, 
to a faster flux decay (a = —0.91). The addition of other available measurements to the RAPTOR 
light curve suggests that another emission component emerged after ~ 10 4 s. We hypothesize that the 
early afterglow emission is powered by extended energy injection or delayed reverse shock emission 
followed by the emergence of forward shock emission. 

Subject headings: gamma rays: bursts - cosmology: observations - shock waves 
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1. INTRODUCTION 

Recent years have brought interesting developments in 
the domain of observations of the early optical emission 
from Gamma- Ray Bursts (GRBs). Several groups now 
have the routine capability to respond to GRB triggers 
in re al time using rapidly slewing robotic instruments 
(e.g. lAkerlof et ab II^TbIoT^EM IBoer. M. 1120041 
Covin p et al. I l2004t iPerez- Ramirez. Park, fc Williams I 
120041 iVestrand et al. 1 l2002j) . Despite much effort in 
this area, so far only a handful of GRBs have been de- 
tected within the first minutes after the onset of the 7-ray 
emission, namely GRBs 990123, 02 1004, 021211, 030418 
041219a, 050319 and 050401 (e.g. fAkerlof et al~l[l999t 
Fox et al. 1120034 iWozniak et al. 1120021 iLi et al. \\200? 
Rvkoff et al. 1120041 IVestrand et al. 1120051 iRvkoff et al. 



2005a,b). Even fewer events have good S/N and cover- 



age. 

The dis covery of the near infra-red transient from GRB 
041219a ijBlake fe Bloom 1120041 iBlake et al 1120051) and 
its p arallel detection in the optical band l|Wren et al. I 
2004) expanded the list of known GRB properties. 
The RAPTO R (Rapid Teles copes for Optical Re- 
sponse; IVestrand et al. I 120021) opti cal light curve of 
GRB 041219a l|Vestrand et a!*1l2005|) overlaps with the 
7-ray emission by an unprecedented ~ 6.4 minutes. 
IVestrand et al. I l|2005f) discovered a qualitatively new 
component of the early optical emissio n from GRBs, and 
prese nted evidence for internal shocks l)Meszaros fc R.eesl 
1999) as the emission mechanism. The presence of the 
new component was established on purely empirical ba- 
sis by its distinct close correlation with strongly time- 
varying 7-ray flux. 

The updated taxonomy for G RB-related optical tran- 
sient (OT) emission proposed bv IVestrand et alTI (2005) 
comprises: (1) Prompt optical emission contempora- 
neous with and consistent with a constant flux ratio 
to 7-rays (the ratio is ~ 1.2 x 10" 5 in GRB 041219a 
ijVestrand et al. 1 12005)): (2) Early afterglow emission 
that may start during the 7-ray emission and lasts for 
several seconds to minutes (uncorrelated with 7-rays 



and typically brighter than the prompt component; e.g. 
GRBs 990123 and 021211); and (3) Late afterglow emis- 
sion that emerges after the fading earl y afterglow and ca n 
persist for hours to many days (e.g. iFox et ai~ll2003aj) . 
The current theoretic al framework offers, cor respond- 
ingly, intern al shocks l| Meszaros fe Reeslll999l). the re- 
verse shock ijSari fe Piran 1 119991 iPan aitescu fe Kumar I 
20041) and the external shoc k ijMeszaros fe Rees 1 119971 
Sari. Piran fc Naravan lll998D phenomena for a possible 
explanation of the observed properties. 

In this letter we present a comprehensive light curve 
of the early optical afterglow emission from GRB 050319 
starting at 35 s after the GRB trigger. 

2. OBSERVATIONS 

On 2005, March 19, 09:31:18.4 UT (trigger time; here- 
after t = 0), the Bur st Alert Telescope (BA T) instrument 
of the Swift satellite l|Gehrels et al. l2004t) detected GRB 
050319, a single-peak event wit h fast raise and exponen - 
tial decay lasting T 90 ~ 10 s l|Krimm et al. ll2005albf) . 
The 15-350 keV fluence, the peak flux and the photon 
index of the time-averaged spectrum were subsequently 
measured to be, respect ively, 8 x 10~ 7 erg cm " 2 , 1.7 ph 
cm- 2 s" 1 and 2.2±0.2 llKrimm et al. Il2005bl) . The on- 
board location ijKrimm etH=bf*"ll2"(~)05aj) was distributed in 
near-real time through the GRB Coordinates Network 
(GCN) at 09:31:36.0 UT, t = 17.6 s. 

Both the RAPTOR-S telescope and the RAPTOR- AB 
array responded to the alert. RAPTOR-S is a fully au- 
tonomous robotic telescope with 0.4-m aperture and typ- 
ical operating focal ratio f/5. It is equipped with a lk 
x lk pixel CCD camera employing a back-illuminated 
Marconi CCD47-10 chip with 13 u pixels. For technica l 
details on RAPTOR- A and B see lVestrand et al. I {2002). 

RAPTOR-S was on target at 09:31:53.7 UT, t = 
35.3 s. The rapid response sequence for RAPTOR- 
S consists of ten 10-second images followed by sixty 
30-second images, a total of ~ 50 minutes of cover- 
age (including the 15-second intervals between expo- 
sures used primarily for readout). A candidate OT at 
a=10: 16:47.9, ^=+43:32:54.5 (J2000) was rapidly iden- 
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Fig. 1. — Examples of single RAPTOR-S exposures of GRB 
050319, the OT discovered by IRvkoff et aT~l I2005al) is circled. 
RAPTOR-S detected the OT with high S/N in the early 10-second 
frames (a, b), and followed its gradual decay down to the magni- 
tude threshold of the 30-second frames (c, d) over the next 40—50 
minutes. The size of the shown area is roughly 4.5' X 4.5'with 
North up and East at the left. 



tified bv IRvkoff et al. I ll2005aft within a hal f hour . The 
OT was later confirmed bv lYoshioka et al~l l)2005(l . and 
the absorption red -shift z = 3.24 was reported by 
iFvnbo et all (2005). Initial analysis of the RAPTOR- 
S images (Fig. showed that the OT was detected at 
high S/N in early exposures and gradually faded be- 
low the magnitude limit. Unfortunately the observing 
conditions at the RAPTOR-S site during response were 
variable and clouds obscured the field of view between 
t = 1480 and 2440 s. 

RAPTOR-B instrument responded slightly faster. Al- 
though none of those images is a detection, including the 
first 10-second frame starting at t = 25.44 s, the corre- 
sponding magnitude limit for OT is of some value (§0}. 
RAPTOR-B and S are separated by about 37 km and 
have independent weather. 

3. PHOTOMETRY 

After standard corrections for bias, dark current and 
flat field responses, all frames underwent a 2x2 pixel bin- 
ning. The binning was applied in order to bring the sam- 
pling of the stellar images to about the critical Nyquist 
value. Additionally, it increased S/N per pixel and made 
the point-spread function (PSF) nearly circular. 

We rejected 29 images taken between t = 1480 and 
2440 s, when transparency was very poor due to pass- 
ing clouds. In even later images the OT detections are 
marginal. We decided, therefore, to form two mean av- 
erages of 5 and 6 frames out of 11 images taken after 
t = 2440 s. A high S/N reference image was prepared 
by mean stacking twenty 30-second frames. All object 
centroids (including OT) were determined using the ref- 



erence image. For the purpose of averaging and photo- 
metric analysis, all frames were resampled to a common 
pixel grid using a bi-cubic spline interpolator and a linear 
coordinate transformation with the ~0.08 pixel accuracy 
(r.m.s.) based on the positions ~ 30 high S/N field stars. 
For image processi ng we used a cu stom Difference Image 
Analysis software ( Wozniak"l l2000j) . 

PSF-weighted photometry within a 4-pixel radius was 
performed assuming a fixed centroid (from reference im- 
age) and without variance weighting. This general tech- 
nique, used in Sloan Digital Sky Survey (Lupton 2005, 
in preparation), hedges against a secondary nonlinear- 
ity between the bright and faint ends of the flux scale. It 
ensures that the much flatter variance profile of the back- 
ground dominated objects cannot propagate the system- 
atic uncertainties from the PSF shape to the photomet- 
ric offsets. We assumed a Gaussian PSF with FWHM 
= 6" (2.44 pixel, binned). The flux scale with about 
3.8% internal consistency was established using 11 high 
S/N stars in the vicinity of the OT. The calibration to 
standard R magnitudes was based on measurements of 
22 USNO-A2.0 stars in the magnitude range R2 = 12.5- 
18.5. Residuals with respect to the best constant mag- 
nitude offset were random over the full flux range (good 
linearity) with r.m.s. scatter of 0.09 mag outside the pho- 
ton noise dominated region. Our unfiltered optical band 
has an effective wavelength close to that of the standard 
R band, but it has a larger width. For lack of the instru- 
mental color information, we assumed that all objects 
have the color of a mean comparison star, i.e. (B — R) = 
1.25 and (R-I) = 0.73, according to USNO-A2.0 cata- 
log. The fact that colors of the early GRB afterglows and 
their temporal evolution remain poorly constrained is a 
source of major uncertainty in transformations of broad 
band photometry (compare §0J. 

The early RAPTOR-B frame was analyzed using the 
same techniques as applied to RAPTOR-S data. The 
actual limit was calculated by performing a fixed centroid 
PSF photometry at numerous random locations near the 
nominal OT position, taking the r.m.s. of the measured 
flux and converting to magnitudes. 

4. RESULTS 

The final RAPTOR photometry of the early optical af- 
terglow of GRB 050319 expressed on the i?-band scale is 
given in Table 2] Fig. |21 plots th e ligh t curve along with 
our model fits. IQuimbv et al. I l)2005fl found an accept- 
able fit to an unfiltered optical light curve from ROTSE- 
Illb telescope using a single power-law model with a — 
-0.59 ±0.05. For the RAPTOR data the best fitting sin- 
gle power-law model has index a = — 0.55±0.02, however 
it yields an unacceptable \ 2 /d.o.f. = 9.20. A visual in- 
spection of the RAPTOR measurements suggests a shal- 
low flux decay at early times and significant steepening 
after ~ 400 s. In fact, the residuals with respect to the 
best fitting single power-law model are systematic and 
indicate a steepening trend. To test that hypothesis we 
fitted a broken power-law model and obtained a reason- 
ably good fit (x 2 /d.o.f. = 2.91) with ai = -0.38 ±0.03, 
a 2 = -0.91 ±0.06 and the break time t hT = 462 ± 55 s. It 
should be noted that instantaneous scale breaking in this 
direction may be hard to explain physically. Neverthe- 
less, we find the model useful for investigating possible 
changes in the light curve slope. 
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TABLE 1 

RAPTOR photometry of GRB 050319. a 
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Fig. 2.— RAPTOR-S optical light curve of GRB 050319 (top) 
and photometric residuals with respect to the reference model (bot- 
tom). Our reference model, i.e. the best fitting single power-law 
(black dotted line), produces systematic residuals. We obtain a 
much better fit with a broken power-law model (red solid line), 
which indicates significant steepening around ~ 400 s. 
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The residuals with respect to the best fitting bro- 
ken power-law model appear flat, however the reduced 
X 2 = 2.91 (81.42/28 d.o.f.) is still formally unaccept- 
able. For some measurements the deviations from the 
best fit model are well in excess of the error estimates, 
in particular there are several strong outliers right near 
the fitted time of the break. There are also a few points 
with fluxes significantly larger than the model prediction 
right before passing clouds covered the field of view. The 
additional photometric scatter could be related to vari- 
able observing conditions and is discussed in more detail 
in m 

To establish the significance of the break, we fitted a 
series of broken power-law models with a range of fixed 
break times. We found that the minimum of the x 2 sur- 
face is not very well constrained and the actual 68% con- 
fidence interval for the break time may be closer to ±100 
s than to the formal parabolic error bar. Our conclusion 
is that the RAPTOR data indicate a significant steepen- 
ing in the flux decay of GRB 050319 within the first hour 
after the 7-ray trigger. However, despite the appearance 
of sharp break in the light curve near ~ 400 s, our data 
are fully consistent with a gradual increase of the slope, 
possibly with additional small scale photometric varia- 
tions. 

5. DISCUSSION 

To test the robustness of the results of §0] we rean- 
alyzed the data using three other photometric tools to 
extract object fluxes: (1) traditional aperture photom- 



a All measurements were obtained with the RAPTOR-S instrument except 
the limit at t = 25 s recorded by RAPTOR-B. Our unfiltered magnitudes 
were transformed to i?-band scale using USNO-A2.0 catalog, and were not 
corrected for extinction (Galactic E(B — V) reddening is only 0.01 mag; 
Schlcgcl, Finkbcinc r fc Davis 1 11998). The last two images are stacks of five 
and six 30-second frames. For those images the effective exposure time ac- 
counts for readout breaks and is shorter than the difference between the end 
and start times (At oxp < t and - i s tart). 



etry; (2) kernel matched difference image photometry 
l Wozniakll2000l): (3) a standard PSF package DoPHOT 
i Schechter. Mateo fc Sahal 11993^1 . and obtained essen- 
tially identical light curves. While the steepening of the 
light curve and most wiggles were always present, the 
precise origin of photometric outliers in Fig. [21 still es- 
capes explanation. Some comparison stars also show the 
wiggles and so me are fu l ly con sistent with the photon 
noise estimate. ILi et alH 1)20031) noticed similar discrep- 
ancies in their light curve of GRB 021211 and suggested 
that color induced systematics could be the cause. The 
intrinsic variability of the OT color or even a stationary 
color difference between the OT and comparison stars 
may generate systematic offsets in photometry. Given 
that the unfiltered spectral band is subject to a red atmo- 
spheric cutoff and weather variations during RAPTOR 
response, it is entirely possible that the outliers come 
from a systematic effect yet to be found. This experi- 
ence underscores the importance of simultaneous color 
measurements of the early GRB afterglows using stan- 
dard filters. 

In Fig. |21 we plot other OT measurements available 
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Fig. 3. — Comparison of the RAPTOR optical light curve of GRB 
050319 with measurements from other instruments. The black 
points are RAPTOR measurements, while the data plotted in color 
are th ose obta ined by other instr uments: ROTSE (Quimbve^al 
|2005l) . UVOT IBovd et al. 120051). Maidanak 1 5 m ^Sharnov et al 
120051). Kiso fe Lulin lm lYoshioka et al. ll2005l) . and ART 14-inch 
ITorii 120051) . The line shows our best fit broken power-law model. 



at this time for comparisons with the RAPTOR light 
curve. The three p oints from ROTSE collaboration 
lIQuimbv et al. l l2005) were shifted by —0.04 mag to re- 
flect the median difference between SDSS r and USNO- 
A2.0 R magnitudes for our com parison stars. Th e V- 
band points from Swift/UVOT ijfBovd et al. ll2005|) are 
plotted 0.7 mag brighter than actual values. All measure- 
ments reported in standard R band were taken at face 
value, since any finer issues with calibration to different 
catalogs should wait until the final revised photometry 
is available. 

Measurements by other experiments agree with the 
shape of the RAPTOR light curve. While more 
sparsely sampled, the V-band light curve observed by 
the Swift/UVOT ijBovd et al. Il2005|) also shows a faster 
flux decay after ~ 400 s. Further, t he me asurement at 
t = 1.27 hours bv lYoshioka et al. I l|2005|) is consistent 
with the extrapolated value predicted by the RAPTOR 
measurements. At times beyond ~ 1.3 hours, the pub- 
lished data from UVOT and other instruments show a 
transition back to a more gradual flux decay rate. The 
two breaks in the flux decay rate are visible in both the 



i?-band and the U-band light curves. The close tracking 
between the light curves in both filters implies a con- 
stant (V — R) color during the first day of the flux evo- 
lution. The zero point fo r the UVOT meas urements is 
still not well established l|Bovd et al. 1 120 05). but taken 
at face value the reported measurements yield a color 
(V — R) = 0.7 for the OT counterpart of GRB 050319. 

While the sample size is still small, one can already 
start to explore the morphology of GRB early afterglow 
light cur ves. The early afterglow behaviors of GRBs 
990123 ijAkerlof et al. I 11999]) and 021211 l)Li et al. I 
2003) were very similar, with both OTs showing a steeper 
initial decline (power-law index a ~ —1.8) and the emer- 
gence of a shallower component (a ~ —0.9) after ~10 
minute s. On the oth er hand, GRBs 0210 04 (jFox et al. I 
I2003b() and 030418 ijRvkoff et al. I [2004) showed shal- 
lower initial decline (or even rise) with a > —0.6 and 
then gradual steepening (to about a < —1.0) on time- 
scales of ~ 10 3 s or longer. The measurements we re- 
ported here for GRB 050319, starting from a\ = —0.38 
and evolving to «2 = —0.91 after ~ 400 s, place its early 
afterglow properties in the latter group. 

In the context of the standard fireball model, the 
shape of the optical afterglow light curve is determined 
by the nature of the interaction between the relativis- 
tic ejecta and the external medium. The relative im- 
portance and timing of the reversed and forward shock 
components, which depend on properties like the den- 
sity profile in the external medium and the strength of 
the magnetic field of the fire-ball, are reflected in the 
rates of flux evolution a nd the break times i n the pre- 
dicte d light curve (e.g. iSari k, Piran 1 119991: iMeszarosI 
l2003 IZhang et al. 112001 . The morphology of the GRB 
990123 and 021211 light curves, with the break to shal- 
lower decay, is usually attributed to the transition from 
the dominance of the reverse shock gene rated emission t o 
forward shock generated emission (e.g. ILi et al. N2003j) . 
For the more gradually decaying early afterglows, the in- 
terpretation is less clear. The gradually declining compo- 
nent could be associated wit h delayed reverse shock emis- 
sion ijVestrand et al. 1 120 05) or an energy injection that 
conti nues well beyond the duration of the initial explo- 
sion (|Fox et al. 12003 rj|) . The emergence of the additional 
component after ~ 10 4 s in GRB 050319 could then be 
understood as the emergence of the forward shock emis- 
sion. The prompt Swift localizations and rapid robotic 
followup are just starting to reveal the richness and com- 
plexity of the GRB afterglow phenomenology and, when 
combined with the models, will help to constrain the ba- 
sic physical parameters of these cataclysmic explosions. 
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